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NONTHERMAL PROPERTIES OF SUPERNOVA REMNANT Gl.9+0.3 
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ABSTRACT 

The properties of the - presumably - youngest Galactic supernova remnant (SNR) Gl.9+0.3 are 
investigated within the framework of nonlinear kinetic theory of cosmic ray acceleration in SNRs. 
The observed angular size and expansion speed as well as the radio and X-ray emission measurements 
are used to determine relevant physical parameters of this SNR. Under the assumption that SNR 
Gl.9+0.3 is the result of a Type la supernova near the Galactic center (at the distance d = 8.5 kpc) 
the nonthermal properties are calculated. In particular, the expected TeV gamma-ray spectral energy 
density is predicted to be as low as e^F^ « 5 x 10~^^ erg cm~^ s~^, strongly dependent [F^ oc d~^^) 
upon the source distance d. 

Subject headings: acceleration of particles — ISM: individual (Gl.9+0.3) — supernova remnants — 
X-rays: individual (Gl.9+0.3) — gamma rays: observations 
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1. INTRODUCTION 

Gl.9+0.3 has been known as a potentially young shell 
type Galacti c supernova remnan t (SNR) of very small 
angular size (jGreen fc Gulllll983). Recently , the interest 
in this SNR was revived bv lRevnolds et al.l (|2008L l2009t ) 
who analyzed the expansion rate of the object and de- 
duced an age ^sn of about 100 yr, which makes it the 
youngest known SNR in the Galaxy. Although the ex- 
pansion rate was derived by a comparison of radio obser- 
vations in 1985 and Chandra observations in 2007, this 
rate has been confirmed very soon thereafter by indepen- 
dent radio observations (jGreen et al.ll2008l : lMurphv et al.l 
[2008). 

According to lRevnolds et al.l ()2008[ ) . the line-free X-ray 
emission has a pure synchrotron origin which clearly in- 
dicates that effective particle acceleration takes place, at 
least for electrons. There are also arguments, like the bi- 
lateral symmetry of the X-ray synchrotron emission sug- 
gesting a roughly uniform ambient magnetic field, that 
favor a type la origin for Gl.9+0.3. Finally, the distance 
estimate d = 8.5 kpc is based on an analysis of the ab- 
sorption toward Gl.9+0.3. 

During the survey of the inner Galaxy by H.E.S.S. 
in very high energy 7-rays, no emission was reporte d 
from the direction to Gl.9+0.3 (jAharonian et aLlbood ). 
Therefore, one can derive an upper limit at the level of 
2% of the Crab flux above 200 GeV. 

For the purpose of a more general study of such 
an unusual object regarding its nonthermal proper- 
ties, it is of interest to describe it by a kinetic the- 
ory of cosmic ray (CR) acceleration in SNRs, cou- 
pled with the gas dynam i cs of the t hermal plasma, as 
given by iBerezhko et al.l (|1996( ) and iBerezhko fc Volfl 
(|1997D . This model assumes spherical symmetry, al- 
though the assumption is later relaxed. Similar models 
on almost the same physical basi s have recently been 
developed by two other groups (|Kang fc JonesI 120061 : 
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IZirakashvili &: Ptuskin I l2009f) , whose calculations very 
well confirm these earlier results. 

The kinetic description allows a corresponding analysis 
of the nonthermal evolution of an SNR at a very early 
phase. This assumes that the plasma physics underly- 
ing especially the temporal dependence of the magnetic 
field amplification process can be extrapolated to such an 
early evolutionary phase, where the dynamical behavior 
of the ejecta plays an essential role. Such a study is re- 
ported here. It is combined with a discussion about the 
influence of the assumption of a smaller distance on the 
TeV 7-ray flux. 

2. MODEL 

Following iRevnolds etHI (|2008D it is assumed that 
Gl.9+0.3 is a Type la supernova (Type la SN) which ex- 
pands into a uniform interstellar medium (ISM). Specif- 
ically, the object is assumed to eject a Chandrasekhar 
mass Moj = 1.4Mq with a total hydrodynamic explosion 
energy i?sN = 10^^ erg. During an initial period, the 
ejecta material has a broad distribution in velocity v. 
The fastest part of these ejecta is de scribed by a pow er 
law dMcj/dv ex. v'^~ ^ with fc = 7 fe.g.. lChevalie3fT98l) . 

The ISM mass density po = l-4TOpiVH, which is usually 
characterized by the hydrogen number density N^i, is 
an important parameter which strongly influences the 
expected SNR dynamics and nonthermal emission; here 
mp denotes the proton mass . 

Following IRevnolds et al.l (j2008') also a distance d = 
8.5 kpc is adopted for the main part of the paper. The 
observed shock size i?s = 2 pc and shock speed Vs — 
14, 000 km s~^ are then used to determine the SNR age 
isN and the ISM number density iVn for the given source 
distance d. 

As reviewed earlier (I Volkl 120041: lBerezhkoll2005l I2008D 
and el aborated most recently in detail in IBerezhko et al.l 
(|2009| ). the key parameters of the theoretical model (pro- 
ton injection flux density, given by a constant injection 
parameter ry <C 1 times the thermal particle flux den- 
sity into the shock, electron-proton ratio below the syn- 
chrotron cooling range, later denoted as K^^, and mag- 
netic field amplification) can be estimated in a semi- 
empirical way from a fit of the theoretical solution to 
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the observed synchrotron emission spectrum if the char- 
acteristics of the initial explosion and the relevant as- 
tronomical parameters are known. For this purpose, the 
nonlinear aspects of the kinetic description are required. 

Of the above-mentioned processes, the magnetic field 
amplification (BelJ, 2004) is the least well understood. 
It is connected with the strong nonlinear excitation of 
magnetic field fluctuations in the shock precur sor by 
the accelerating energetic particles ( McKcnzie fc VolU 
fl98llLucek fc Belll[2000l: iBell fc Lucekll2001l: IBelil200ir 
Therefore, it is assumed here that these fluctuations 
lead to Bohm diffusion of the energetic particles in 
this amplified field. Such a bootstrap mechanism 
can be approximately ju stified by the resu lts of re- 
cent particle simulations ([Reville et al.l [2008[ ). In ad- 
dition, the amplification process is strongly dissipa- 
tive, as shown by hydromagnetic and k i netic simula- 
tions ('Beli"200l: 'Zirakashvili e t al.l [200I: iReville et al.l 
2008; Nicmicc et al. 2008; Riaue lme fc Spitkovskvl 120091 : 
Ohira et al.ll2009l ). Therefore, it must be accompanied 
by strong gas heating within the precursor region r > Rg 
due to wave dissipation, adopted in the present model 
in the form (def,/dt)r \ ,ss = —a'»c\dPr.l d r with an = 1 
(|Berezhko et al.ni99alBerezhko fc V61klll997D . where eg 
is the gas thermal energy density, Pc denotes the ener- 
getic particle pressure, and where, in a second bootstrap 
mechanism, ca = i?(47r/9)^^/^ is the Alfven velocity in 
the amplified field B (see below). The value an = 1 cor- 
responds to the assumption that the Alfven wave field ex- 
cited within the precursor reaches amplitudes which are 
very much smaller than the maximal amplitudes which 
could be reached if the whole work —cpJdPd dr done by 
CRs went into wave excitation. In such a case, this work 
goes almost completely into gas heating due to the wave 
damping. The gas thermal pressure just ahead of sub- 
shock is in this case consi derably larger th an the pres- 
sure of the magnetic field ([Berezhkol I2008D . Therefore, 
the subshock can be treated approximately as a pure 
gas shock. This approximation will later be re-examined 
through the approximate inclusion of the amplified field 
and its associated turbulent gas motions in the subshock 
dynamics. 

The magnitude of field amplification in all young SNRs 
is such that a non-negligible fraction of the shock en- 
ergy PnK^ is co nverted into magnetic field energy (e.g. 
iBerezhkol I2008D . In fact, a time dependent, amplified 
upstream magnetic field strength 

B^(t) = BoitsN)[Vs{t)/VsitsN)]' (1) 
is used here, where the theory parameter i3o(isN) is the 
rms field strength at the present epoch isN and is es- 
timated by a comparison of the theoretically calculated 
synchrotron spectrum with the observed one (see below). 
Such a form of the time dependence of the amplified mag- 
netic field with (5 « 1 is consistent with the interior field 
strengths estimated from observational resu lts for a num- 
ber of young SNRs (e.g. iVolk et al.l[2005h . (Since iBelll 

(I2OOI even estimated a dependence Bo{t) oc Vs^^(t) for 
the field amplification due to the nonresonant stream- 
ing instability alone, also the case of 5 = 3/2 will be 
examined here. The radial dependence of the rms. mag- 
netic field strength B{r,t) in the shock precursor is then 
modeled by B{r,t) — Bo{t)p{r,t)/ pn, where po is the 
far upstream (interstellar) density. In the overall con- 



servation relations for momentum and energy of the sys- 
tem, the magnetic field strength in the upstream ISM is 
BisM < Bq. Thus, Bo/BisM is the field amplification 
factor by the accelerating energetic particles alone. An 
anal ysis per formed for a number of young SNRs shows 
(e.g.. lBerezhko. 2008) that the field strength Bq, required 
to fit the observed synchrotron spectrum, is well within 
the range expected from t heoretical estimates fe.g. jBell 
[2OOI: IPeiietier et al.l[2006h . 

The observed X-ray morphology of SNR Gl. 9-1-0. 3 
agrees with the theoretical expectations regarding the 
morphology of ion injection and the corresponding mor- 
phol ogy of magnetic field amplification for a Type la 
SN (jVolk et al.l I2003D . It is therefore consistent with 
a correction for the spherically symmetric solution by 
a renormalization factor f^c ~ 0.2 of the energy den- 
sity of nuclear partic les, like in the case of SN 1006 
(|Berezhko et al.ir2009( ). In this case, the electron-proton 
ratio Kep, calculated under the assumption of spherical 
symmetry, should be increased by a factor l//ic- 

3. RESULTS 

The calculated evolution of the gas dynamical vari- 
ables of Gl. 9-1-0. 3 is shown in Figure [TJ The observed 
shock radius Rg and shock speed Vg are fitted at the age 
tsN = 80 yr and for an ISM hydrogen number density 
A'h — 0.018 cm~^ (Figure [1^). Since the remnant is 
in the free expansion phase, it is approximately consis- 
tent with an an alytical self-similar solution Rg oc t^/^ 
(jChevalieil Il982[ ). Note that an explosion model with 
an exponential ejecta velocity profile gives slightly larger 
values of the age, tsN = 100 yr, and of the ISM density, 
corresponding to a hydrogen number density of about 



0.03 cm" 



for the same assumed distance of 



d = 8.5 kpc (jRevnolds et al.ll2008D . 

The calculated radius Rc of the contact discontinuity 
(CD) and the CD speed Vc are also shown in Figure[lja). 
One can see that the ratio Rc/Rs is rather small. At the 
current epoch Rc/Rs ~ 0.9. 

To obtain a good fit for the observed synchrotron 
spectrum (see below) first of all a proton injection rate 
rj = 10~^ is required. It leads to a nonlinear modification 
of the shock which, at the current age of t = 80 yr, has 
a total compression ratio a w 4.6 and a subshock com- 
pression ratio CTs ~ 3.6 (Figure [TJb)). In addition, an 
electron-proton ratio Kep « 5 x 10~ and an upstream 
amplified magnetic field strength at the current epoch 
i?o(isN) ~ 100 pG are required. This implies a down- 
stream magnetic field strength of B^ ~ 460 pG. 

All the above quantities are almost the same for the 
two different time dependences of the magnetic field 
B(){t) within the evolutionary period t = igN i 100 yr. 
The only difference between these cases is a slightly more 
rapid increase of the shock modification (characterized 
by the shock compression ratio a{t)) in the case 5 = 3/2 
compared with the case 6 = 1 due to the increase of the 
Alfvenic Mach number Ma oc Vs/Bq oc V^^""^ oc i3(<5-i)/7_ 

To explain these results it is noted that, as in other, 
similar cases of strong, modified shocks, the existing 
measurements permit an estimate of the three parame- 
ters of the theoretical model. This takes into account 
the following influence of the parameters on the syn- 
chrotron spectrum. (1) Since the high-energy electrons 
undergo strong synchrotron losses and, since during the 
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Fig. 1. — Shock (CD, (dashed Unes)) radii _Rs (iJc) and shock 
(CD, (dashed Unes)) speeds Vs (Vc) in units of _Ro = 1 pc and 
Vb = 10'* km s~^ (a), total shock a and thermal subshock 
compression ratios (b), and total energy contents of accelerated 
CRs, Ec (c), as functions of time in years. The red colored curves 
imply an amplified magnetic field strength B cc Vs(t), whereas 

the blue colored curves correspond to i? oc Vs^'^(t), cf. Equation 
(1); the difference between these two cases is not distinguishable 
in Figure [TJa). The dotted vertical line marks the current epoch, 
iSN = 80 yr. The observed mean size (open circle) and speed 

i filled circle) of the sh ock, as determined by X-ray measurements 
Reynolds etalll2008h . are shown as well. 

acceleration process they dominate the nonthermal elec- 
tron pressure P°, the flux of X-ray synchrotron emis- 
sion which they produce vSi, is approximately propor- 
tional to the total energy flux of nonthermal electrons 
Fe oc P^VsRl oc KcpPcVsR^, and is only weakly sensitive 
to the magnetic field strength _Bo- The spectrum of high- 
energy protons N{p), which give the main contribution 
to the total CR pressure Pc ^ pV^, is only weakly sen- 
sitive to the injection parameter 77 in the case of a mod- 
ified shock. Therefore, the fit of the observed X-ray flux 
mainly determines the value K^p of the electron:proton 
ratio. The spectrum of accelerated electrons is calcu- 
lated in absolute numbers. It is expressed in terms of 



Kcp, due to the presumably dominant dynamical role of 
protons. (2) Values a > 0.5 of the radio spectral index 
a — —dlnSi^/dlnv, as observed in young SNRs, require 
a modified shock with CTs < 4 < cr (This also implies a 
curved electron spectrum that hardens toward high fre- 
quencies). The value of a is mainly determined by the 
subshock compression ratio Us, which in turn is deter- 
mined mainly by the proton injection rate 77. Therefore 
the fit of the measured spectral shape of the radio syn- 
chrotron emission gives mainly the required value of the 
injection parameter rj. (3) Finally, since the radio emis- 
sion flux value Si^ oc Kcpi'^°' B'^^^ is strongly dependent 
upon the magnetic field strength, its value is derived 
from the fit to the observed amplitude of the radio spec- 
trum. Thus, three measured characteristics of the syn- 
chrotron spectrum — X-ray fiux, shape, and amplitude 
of the radio emission — make it possible to obtain an 
estimate of the three relevant "theory parameters" B^, 
Kcp and rj even though this is not a simple three-step pro- 
cedure but an itera tive procedure, minim izing the com- 
bined x^-value fsee iBerezhko et al1l2009l for details). 

Note that the magnetic field value B^ can also be es- 
timated by another, independent method, namely from 
a fit of the observed spatial fine structure of the X-ray 
emission. In all the cases where such measurements 
exist, both meth ods give consistent values of B^ (e.g., 
IVolk et all 120051 ). Unfortunately, the fine structure of 
the X-ray emission is not determined yet for Gl. 9-1-0. 3. 
Therefore, this consistency check cannot be made at 
present. 

The uncertainties of the estimated values of rj, i?d, and 
Kcp depend upon the quality of the measurements of 
the synchrotron spectrum and can be rather small, as 
it was recently demo nstrated for the case of SN 1006 
(|Berezhko et al.l[200 9^. In the case under consideration 
it is about 20% for and Kop and 30% for 77. 

It should also be noted that due to the very low age 
of the SNR and the low ISM density the expected ther- 
mal X-ray emission is far below the observed X-ray fiux, 
which is therefore completely dominated by the nonther- 
mal component. 

The required proton injection rate rj = 10~^ is consid- 
erably higher than the critical value ry* w 10"'*, which 
separates a nonlinearly modified shock with 77 > 77, from 
an unmodified st ate, resulting fr om very low injection 
rates 77 < 7/* fsee IBerezhko fc Elli son 1999, Eq.38). The 
relatively weak shock modification is the result of the 
very large magnetic field Bq, that leads to strong gas 
heating within the precursor region r > Rg. 

In order to check the sensitivity of the results to the 
adopted value of the parameter an, calculations with 
0.02 < Oh < 1 were performed. It turned out that even 
for an — 0.02 the shock properties are not very far from 
the case an = 1: the required far upstream magnetic 
field value is Bq = 125 /iG, and the shock compression 
ratios are a = 6.3 and ag — 3.99. 

It is noted here that the results of IVladimirov et al.l 
(|2008l ) suggest a stronger sensitivity of the shock proper- 
ties to the value of an . The reason is that in the consider- 
ations of these authors the parameter an determines not 
only the gas heating through the term —anCAdPc/dr, 
but also the upstream magnetic field generation through 
a complementary term (an — \)cxdPc/ dr. Such an ap- 
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proach implies that magnetic field amplification takes 
place only due to resonant excitation of Alfven waves 
and, in addition, according to such a quasi-linear ex- 
pression. In contrast, the present model allo ws that the 
magnetic field is also amplified nonresonantly ('Bel]| |20d5 
iPelletier et al.l l2006t ) . Secondly, the determination of 
the amplified field i?o(^SN) in the present semi-empirical 
model uses the synchrotron observations of the source. 
For both these reasons, the magnetic field amplification 
and the gas heating are not connected by a simple rela- 
tion. 

3.1. Subshock dynamics with amplified B-field 

([1981 



Following 


McKenzie & Volkl 


iVladimirov et al. 


( 20081) ■ an I 



and 
connection 

w = B(47rp)~-'^/^ between the magnetic field fluctuation 
vector B and the fiuctuation vector w of the mass 
velocity is assumed. This implies an approximately 
incompressible plasma turbulence with a locally ho- 
mogeneous mass density p, where the total (mean 
square magnetic field plus plasma turbulent) pressure 
is given by i-turb = (B)^/(87r), and the total (mean 
square magnetic field plus plasma turbulent) energy 
flux density equals Jlurb = SuPturbi here u denotes the 
shock-normal mean mass velocity in the shock frame. 
These normal components of the momentum and energy 
flux densities, immediately upstream and downstream 
of the subshock were included in the Rankine-Hugoniot 
conditions for the subshock. The latter is approximated 
as a locally plane, normal shock wave. (In the precursor 
region, (B)^ = is the mean square strength of the 
amplified magnetic field B, introduced in the previous 
section and assumed here to be isotropically distributed 
with Gaussian statistics; the field strength downstream 
of the subshock is taken to be B2 — <JsBi, where Bi is 
the field strength upstream of the subshock.) Note that 
the above expression for the total turbulent energy flux 
density Jturb in the downstream region differs somewha t 
from the expression used by ICaprioli et al.l ()2008[ ). 
Since the consideration of these authors is based on the 
transmission and reflection of small-amplitude Alfven 
waves at purely parallel subshock, it is believed here 
that this linear treatment is not applicable to the actual 
case of a strongly perturbe d and amplified niagneti c field 
(see also the arguments of IVladimirov et all (|2008l )). 

It is clear that the approximations introduced above do 
not exactly describe the true physical situation thatalso 
contains th e nonreson antly unstable modes of the Bell 
instability (jBelll 1200^ , because these transverse modes 
will in their nonlinear evolution also develop compress- 
ible elements (BcU 2004; Zirakashvfli et al. 200^^ In 

addit ion, also the acoustic modes (Dorfi 19841: iDrurvl 
11984 lBerezhkolll986t iMalkov fc Dia^ond.2G06i rwifl con- 
tribute. Their influence on the strength of the subshock 
remains to be evaluated. However it is believed that the 
present description gives a roughly correct estimate of the 
subshock effects of at least the incompressible part of the 
fluctuation fields produced by the accelerating particles. 

In the adopted approximation a — 4.55, Ug — 3.6 for 
Q^H = 1, very close to the previous case, where the tur- 
bulent momentum and energy fiuxes were ignored in the 
subshock conservation relations. The effect is somewhat 
larger for the smallest value of the parameter an = 0.02: 
it leads to a decrease of a from a — Q.i to a — b.'d 



and to a decrease of the magnetic field strength from 
i?o — 125 jiG to So = 120 /iG, which is not a large effect 
either. 

It should also be noted that the assumption of consid- 
erable gas heating due to wave dissipation, correspond- 
ing to QfH = 0.5 — 1, is consistent with the numerical 
modeling of the nonres onant wave excitation (lBel]|l2004l: 
IZirakashvili et al.ll2Q08[ ). Such dissipation should operate 
in a similar way for the resonant Alfven mode instability. 
It is therefore concluded that insignificant gas heating, 
which occurs for an ^ 1 within the present formalism, 
is an unrealistic assumption. 

3.2. Charged particle and ^-ray spectra 

With the renormalization /ic = 0.2, the nuclear CRs 
inside Gl.9-1-0.3 SNR contain (Figure ffljc)) 

w 0.0025£;sN ~ 3 X 10"*^ erg. (2) 
The volume-integrated (or overall) CR spectrum 



7V(p,0 = 167rV / drr^f{r,p,t) 



(3) 



has, for the case of protons, almost a pure power- law form 
N (X p~^ over a wide momentum range from O.lmpC up 
to the cutoff momentum Pmax ~ 3 x lO^mpC (Figure [2]). 
This value Pmax oc RsVsBq is limited mainly by the fi- 
nite size and speed of the shock, its deceleration and the 
adiabatic cooli ng effect in the downstream region (see 
lBerezhkolll996l for details). As pointed out above, par- 
ticle diffusion is approximated by Bohm diffusion in the 
amplified magnetic field B, cf. Equation (1). It is impor- 
tant to note that the calculated value of Pmax is therefore 
an upper limit, because it is assumed that up to the cut- 
off all particles "see" the amplified field everywhere. 

Consequently, Gl.9-1-0.3 represents the youngest SNR 
where the accelerated proton spectrum extends up the 
so-called knee energy. Such a maximum proton energy 
appears indeed required to desc ribe the overall CR spec - 
trum for energies up to 10^"^ eV (jBerezhko fc Volk|[2007l) . 

The shape of the overall electron spectrum (p) devi- 
ates from that of the proton spectrum N{p) at high mo- 
menta p > p\ ^ lO'^mpC on account of the synchrotron 
losses during the electron residence time in the down- 
stream region (Figure [2]) . Within the momentum range 
p\ < p < pjjiaxj the electron spectrum is considerably 
steeper, Nf. oc p~^, due to synchrotron losses taking 
place in the downstream region after the acceleration 
at the shock front. The maximum electron momentum 
Pmax ~ lO^TOpC corrcspouds closely to the result obtained 
by equating the synchrotron loss time and the accelera- 
tion time. 

Figure [3] illustrates the consistency of the synchrotron 
spectrum, calculated for the above-mentioned best set of 
parameters with the observed spatially integrated spec- 
tra. 

As mentioned above, values a > 0.5 of the radio spec- 
tral index a = —dlnS^/dlnu, as observed in young 
SNRs, require a curved electron spectrum that hardens 
toward higher energies, as predicted by nonlinear shock 
acceleration theory. To hav e a = 0.62 in the ra dio range, 
as observed for Gl.9+0.3 ([Green et al.l l2008\ requires 
efficient CR acceleration with a proton injection rate 
rj — 10~^ which leads to the required shock modification, 
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Fig. 2. — Spatially integrated CR spectrum as function of particle 
momentum. Solid and dashed lines correspond to protons and 
electrons, respectively. 

and also leads to the high magnetic field value above. As 
it is clear from Figure [3l a good fit of the observed X-ray 
energy flux can only be achieved due to the softening of 
the synchrotron spectrum for ^>10^^ Hz, which is due 
to the strong synchrotron losses of electrons with mo- 
menta p > pi ~ 700mpC. Using the known dependence 
pi oc B^^t''^, a.t t = isN one can immediately estimate 
the required value of the interior magnetic field value 
Bd ~ 500 /iG, consistent with the above determination 
from the radio spectrum. The calculated hard power-law 
X-ray spectrum continues almost up to 50 keV, making 
this source in principle attractive to be observed with 
Suzaku and INTEGRAL and the future Astro-H X-ray 
instrument. 

A very important question for every young SNR is 
whether the existing data indeed unavoidably require ef- 
ficient proton acceleration accompanied by strong mag- 
netic field amplification. In order to explore the alter- 
native possibility, Figure [3] presents in the dotted curve 
a synchrotron spectrum which corresponds to a hypo- 
thetical leptonic scenario with a proton injection rate 
so small [r] <^ lO^*^) that the accelerated nuclear CRs 
do not produce any significant shock modification and 
therefore also no magnetic field amplification. This cor- 
responds to the test particle limit, when the distribution 
function of shock accelerated electrons has the form 

/o = Ap^^ exp(-p/pniax), (4) 

where the amplitude A and the value of the cutoff mo- 
mentum Prnax are determined by the fit to the observed 
synchrotron spectrum for a given interior magnetic field 
value Bd- Since magnetic field amplification is not 
expected in this case, the downstream magnetic field 
cannot be larger than the MHD-compressed ISM field 
-BisM ~ 5 /^G. The maximal possible downstream field 
-Bd ~ 20 is adopted which corresponds to the mini- 
mal number of accelerating electrons, and therefore the 
7-ray emission produced by these electrons is also min- 
imal. The synchrotron spectrum for the leptonic test 
particle scenario in Figure [3] corresponds to the maxi- 
mal electron energy ej^jj^x = -Pmax"^ — 6 TeV, determined 
by the Chandra observation. There are two differences 
in the synchrotron spectra, corresponding to these two 
scenarios. The high-injection scenario leads to a soft ra- 



Chandro 




a 10 12 14 16 IS 20 

log v, Hz 



Fig. 3. — Spatially integrated synchrotron SED as a function 
of frequency. The dotted line corresponds to the test particle 
limit (leptonic scenario). Fluxes of the X-ray emission observed 
by Chandra (Reynolds ct al. 2008) and the radio emission com- 
piled by Green ct al. ( 2008) are also shown. The solid lines are 
fitted to the mos t rece nt VL A rad io data shown by the open cir- 
cles, reported bv I Green et aTl (120081 ) . The line colors have the same 
meaning as in Figure [T] 

dio spectrum Si, oc v^"' with power law index a — 0.62, 
whereas in the test particle case a = 0.5. On the other 
hand, the two spectra behave essentially differently at 
X-ray frequencies J^^IO"'^^. This demonstrates that only 
in the high-injection case with its high, amplified mag- 
netic field value Bd « 460 /iG the spectrum S^iy) has a 
smooth cutoff, consistent with the observations (see Fig- 
ure [3]). In the test particle case the spectrum S^iy) has 
too sharp a cutoff to be consistent with the observations. 

If Gl. 9-1-0. 3 was indeed a Type la SN, then the explo- 
sion parameters -Esn, -^cj and k are known. Assuming 
the value d = 8.5 kpc for the source distance and us- 
ing the age tsN and the ambient gas number density N^i 
from fits to the observed astronomical parameters size 
and expansion rate, and also using the "theory parame- 
ter" values ?7, K^p and -Bq, estimated from the fit to the 
synchrotron spectrum, one can predict the 7-ray fiux for 
the assumed source distance d. 

In Figure IH the calculated 7-ray spectral energy dis- 
tributions (SEDs) due to Tr'^-decay and inverse Compton 
(IG) collisions are presented for the source distance d — 
8.5 kpc together with the sensitivities of the Fermi and 
H.E.S.S. instruments. For the modified shock, consistent 
with the observed synchrotron emission, the expected to- 
tal TeV 7-ray SED is e^F^ « 3x lO'^ eV cm'^ s'^ Such 
a fiux is too low for an H.E.S.S. detection in ^ 50 hr by a 
factor of the order of 30. The Tr'^-decay flux is only about 
6% of the IC 7-ray flux as a result of the low gas density . 
Since according to recent estimates ([Porter et al.ll2006[l . 
the Galactic interstellar optical and infrared radiation 
fields in the inner Galaxy are considerably higher than 
previously thought, for this region the calculation of the 
IG fiux was performed on the basis of those estimates. 
The higher radiation field leads to an increase of the IG 
gamma-ray flux at energies e-y < 1 TeV by an order of 
magnitude compared with a standard inte rstellar radi- 
ation fleld in the solar neighborhood (e.g. iDrurv et al.l 
119941 ). Nevertheless, as indicated above, the expected 
TeV emission flux is still far below the H.E.S.S. sensitiv- 
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ity. 

As can be seen from Figure lU the TeV 7-ray flux ex- 
pected in the unmodified leptonic scenario considerably 
exceeds the H.E.S.S. sensitivity, corresponding to ^ 50 hr 
of observation time. Since the region of the Galactic cen- 
ter was already explored by H.E.S.S. for times of more 
than 100 hr without detection of Gl.9+0.3, this purely 
leptonic test particle scenario should be rejecte d as in all 
similar cases of Type la SNe (jVolk et al.ll2008[ ). 

3.3. Dependence on the assumed source distance 

It is, however, to be noted that the expected 7-ray flux 
is very sensitive to the assumed source distance d. There- 
fore, SNR Gl. 9-1-0. 3 could be a potential 7-ray source if 
the actual distance was lower than 8.5 kpc. Qualita- 
tively, the dependence of the expected 7-ray flux on dis- 
tance can be understood if one takes into account that 
the 7r°-decay 7-ray energy flux e^F^ oc MswGc/d'^ is pro- 
portional to the mass of gas, Mgw = 47ri?gPo/3, swept up 
by the SN shock, and to the energy density Cc of the CRs 
producing 7-rays of given energy. Since for high accel- 
eration efficiency Cc is proportional to the shock kinetic 
energy density poV^, one can write 

e-yF^ N^V^R^/d^. (5) 

For fixed explosion energy Esn, the distance d and ISM 
density A'h are connected by the relation 

A^H oc (6) 

because in the free expansion phase the SNR radius 

i?s oc is determ ined by the expression Rs oc N^^^'^t^^'^ 
(|Chevalieijll982l ). where the SNR age t oc Rs/V^ is fixed if 
the angular size and angular expansion speed are known 
as in our case of Gl.9+0.3. 

Taking also into account that for a fixed angular ex- 
pansion rate of the object Vg cc d, one obtains 

e-yF-y(xd-^\ (7) 

According to this relation a mere 30% reduction of the 
source distance leads to an increase of the expected 7- 
ray flux by a factor of more than 10. This is illustrated 
in Figure m where also 7-ray spectra are presented that 
were calculated for the distance value d = 5.6 kpc. In 
this case the shock velocity and size could be fitted at 
the same age tgN = 80 yr and for an ISM hydrogen num- 
ber density A'h — 0.2 cm~'^. A similar fit for the radio 
and X-ray data as in Figure |4] could be achieved with 
an electron-proton ratio K^p = 8 x 10~* and a down- 
stream magnetic field strength = 670 /iG. It is clear 
that Gl. 9-1-0.3 could be visible in TeV 7-rays by future 
instruments like the Cherenkov Telescope Array (CTA), 
if the actual distance was not larger than d = 5.6 kpc. 

The theoretical, spatially integrated radio synchrotron 
flux slowly increases with time, as can be seen in Figure 
m essentially due to the rapidly increasing total number 
of accelerated electrons in the increasing SNR volume 
cx Rl 

The X-ray synchrotron flux is expected to be nearly 
constant in time (Figure [5]). This is mainly due to the 
strong synchrotron cooling of the highest energy elec- 
trons which produce the X-ray synchrotron emission. 

The TeV 7-ray flux is expected to increase with time 
as well (Figure [5]), mainly due to the increase of overall 
number of GRs with energy above 10 TeV. 
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Fig. 4. — Integral Tr'^-decay (solid lines) and IC (dashed lines) 7- 
ray energy flux densities (SEDs) for the present epoch, as a function 
of 7-ray energy, for the two different source distances d = 8.5 kpc 
(thick lines) and d = 5.6 kpc (thin lines). The calculations are given 
for the case 5 = 1. SEDs calculated for the case <5 = 3/2 coincide 
with the presented SEDs within 10% accuracy. The thick dashed 
line represents the IC 7-ray energy flux, calculated for the recently 
re-estimated interstellar optical/infrared radiation background in 
the central region of the Galaxy ( Porter et al! 200l). The dotted 
line corresponds to the test particle limit (leptonic scenario). For 
comparison, the sensitivities of Fermi (for a 5cr detection in 1 yr 
of sky survey exposure with a background representative of the 
diffuse background near the galactic plane; lAtwood et al.l {2009) ) 
and H.E.S.S. (for a 5cr detection of the Crab Nebula with powcr- 
law d i fferen tial photon index 2.6 in 50 hr at a zenith angle of 20°; 
[Fuiiia 1I2OO5M are shown. 
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Fig. 5. — Time dependence (in years) of the fluxes of the ra- 
dio synchrotron emission at frequency v = 1.4 GHz (dash-dotted 
lines), synchrotron X-ray emission with energy = 6 keV (dashed 
lines), and TeV energy 7-ray emission (solid lines). The fluxes are 
normaliz ed to thei r values _Fo at the current epoch. Available data 
in radio I) Green et al. 2008) are shown as well. The line colors have 
the same meaning as in Figure [T] 

3.4. Renormalization of the nuclear particle spectra 

If the ambient interstellar magnetic field would be 
completely disordered on spatial scales smaller than the 
shock size, then efficient CR injection/acceleration would 
be expected across the whole shock surface. In such a 
case /re ~ 1 and the expected gamma-ray flux would 
be higher by a factor of l//ro ~ 5, whereas all the 
emission produced by GR electrons would remain the 
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same since the normahzation of the electron spectrum 
was done based on the observations. However, it is be- 
lieved that the actual situation is opposite: the bilateral 
symmetry of the X-ray synchrotron emission suggests a 
roughly uniform ambient magnetic field on a parsec scale 
and therefore Ac « 0.2 like in SN 1006. 

4. CONCLUSIONS 

The existing data for Gl.9+0.3, when analyzed within 
the framework of the nonlinear kinetic theory of CR pro- 
duction in SNRs described above, are consistent with a 
type la explosion in a rarefied medium at a distance of 
d = 8.5 kpc, whose nuclear CR spectrum reaches the 
energy of the "knee" in the observed Galactic CR spec- 
trum at the present epoch. This conclusion also concerns 
the derived strong magnetic field amplification. A test 
particle, purely leptonic gamma-ray scenario, is inconsis- 
tent with existing TeV gamma-ray observations with the 
H.E.S.S. telescope array. However, the data set is not 
complete enough to unequivocally determine the value 
of the source distance. Since the expected 7-ray flux 
is very sensitive to the distance d, oc d^^^, a detec- 
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it may be, would yield the distance. However, in the 
case when Gl. 9-1-0. 3 is located near the Galactic center 
{d = 8.5 kpc), the expected TeV 7-ray energy flux is so 
low, e^F^ « 5 X 10~^^ erg cm~^ s~^, that it is not de- 
tectable with present instruments. It is clear that for a 
distance that was not larger than d = 5.6 kpc Gl. 9-1-0. 3 
could be visible with future instruments like CTA which 
can be assum ed to have a se nsitivity of ^ 1 mcrab at a 
few 100 GeV (jBernlohijl2009l ). 
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